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I. INTRODUCTION 

Transitions with a change of Euclidean topology is a 
subject of wide physical interest. One interesting exam- 
ple is the phase transition connected with a nucleation 
of a black hole in a thermal bath. Consider a thermal 
field with temperature T in a flat spacetime. One can 
use Euclidean fields on a spacetime with the topology 
R ^ 1 x S* 1 to describe a canonical ensemble for such a 
field. The size of the compact dimension S 1 is ft = l/T. 
A nucleation of a black hole changes the Euclidean topol- 
ogy from R D ~ 1 x 5 1 to S D ~ 2 x R 2 . The corresponding 
Euclidean space after the black hole nucleation is the 
Gibbons-Hawking instanton 

Another important example of a similar phenomenon 
is the so-called merger phase transition which occurs in 
models with large extra dimensions when a black hole is 
localized in a spacetime which has additional k compact 
dimensions (D = 4 + k) (a caged black hole, for reviews 
see @, H, 0, M, El)- In the absence of the black hole such 
a spacetime has the topology R 4 x T k . Kol argued Q 
that the black hole-black string phase transition includes 
a local topology change of the corresponding Euclidean 
manifold so that the singular geometry is a cone over 
S D ~ 3 x S 2 8]. This topology change is similar to the 
conifold transition Q. In the black hole phase the S D ~ 3 
is contractible while in the black string phase the S 2 is 
contractible. In order to achieve this topology change 
one has to pass a configuration which is singular at the 
tip of the cone. The "double-cone" over S D ~ 3 x S 2 is 
given by 

ds 2 = dp 2 + ^-^ [d X 2 + cos(x) 2 dt 2 + (D- 4) d02,_ 3 ] . 

__ (1) 

For more details see [2, LL0| • 
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Kol [Tl| proposed that there exists a relation be- 
tween merger transitions and Choptuik's critical collapse 
[i~2l [l3j . This correspondence can be achieved by per- 
forming two analytic continuations. The physics of the 
critical collapse and merger transitions have some com- 
mon features like a singular critical solution which turns 
out to be an attractor and a self-similar solution in the 
neighborhood of the singular point. A better understand- 
ing of one of the systems may shed light on the other. 

It is interesting that there also exists a close similarity 
between the properties of merger transitions and a toy 
model proposed some time ago for study transitions dur- 
ing which the Euclidean topology is changed [3, EH, EH • 
This model consists of an A-dimensional static bulk black 
hole and a D-dimensional brane (D < N) interacting 
with this black hole. The brane is assumed to be a test 
brane and infinitely thin. The former assumption means 
that one neglects the effects connected with the gravi- 
tational field of the brane, while the latter one implies 
that the effects of the brane thickness are neglected and 
its worldsheet is a minimal surface which provides an ex- 
ternum of the Dirac-Nambu-Goto (DNG) action. It is 
assumed that the brane is static and axisymmetric, so 
that the induced geometry on the brane possesses the 
0(D — 1) group of isometry. It is also assumed that 
far from the black hole the brane surface is parallel to 
the equatorial plane of the bulk black hole. For such a 
brane there exists two qualitatively different configura- 
tions: One, which is called subcritical, is a brane which 
does not intersect the black hole event horizon, and the 
other, supercritical, is a brane crossing the horizon. In the 
latter configuration the induced geometry on the world 
sheet of the brane is the geometry of a D-dimensional 
black hole, which is called a brane black hole, or briefly 
BBH. Such a black hole is absent for a subcritical con- 
figuration. Thus by changing the position of the brane 
at infinity (asymptotic data), one generates a transition 
between BBH and no-BBH phases (see Fig. [TJ. 

If this change is done adiabatically, then one deals with 
a one parameter set of quasistatic solutions. After Wick's 
rotation of time one gets a one parameter set of Euclidean 
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FIG. 1: Three possible types of configurations - the subcritical 
embedding (left) where the brane does not touch the black 
hole horizon, the supercritical (right) where we have an in- 
duced black hole on the brane (BBH), and the critical (center) 
which is singular at its tip. 



induced metrics, with a change of the Euclidean topol- 
ogy of the induced metric at some critical value of the 
asymptotic data, which plays the role of an order param- 
eter. The Euclidean topology changes from S 1 x BP -1 
for the subcritical configuration to R 2 x S D ~ 2 for the su- 
percritical. It was demonstrated [3, EH EH] that when 
the effects of the brane stiffness are neglected the rela- 
tion between the asymptotic data and the mass of the 
induced BBH for the transition between sub- and su- 
percritical configurations is universal, that is it does not 
depend on the bulk black hole characteristics. Moreover, 
there is no mass gap for a creation of a BBH, so that 
the corresponding phase transition is of the second order. 
The near-critical solutions possess discrete (for D < 6) 
or continuous (for D > 6) self-similarity, which makes 
this transition formally similar both to the merger tran- 
sitions and the near-critical collapse discovered by Chop- 
tuik 1131 . (For a general review of the critical collapse see 
e.g. [13| |. See also a discussion [TtJ of the critical collapse 
in a higher dimensional spacetime.). These properties 
and close similarity of near-critical solutions for both the 
BBH model and merger transitions make it interesting 
to consider in the framework of the BBH model some 
general problems which exist for this class of models. 

Before discussing these problems we mention that the 
universality of the near-critical behavior in the BBH 
model is a consequence of the following fact: only near- 
horizon properties of these solutions are important . In 
this near-horizon domain there is no dimensionful param- 
eter which determines the behavior of the system. As a 
result of this the system has scaling properties and the re- 
lated phase transition is of the second order. This is why 
the model captures many universal features of various 
physical systems (l8| . One important example of such a 
system provides a holographic description of the meson 
melting phase transition of matter in the fundamental 
representation [H, HI US HI HI- The configuration 
consists of N c color Dp-branes and Nf flavor Dq-branes 
when p < q. The addition of the flavor Dq-branes is 
dual to the addition of matter in the fundamental repre- 
sentation in the gauge theory. In the limit N c ^> Nf the 
Dp-branes are described by a p-brane supergravity action 



("black Dp-branes") while the Dq-branes are described 
by the DNG action. In other words, we can say that 
there are Dq-branes in the background of Dp-branes. In 
this system the point, where the brane touches the black 
hole horizon corresponds to a certain temperature where 
the mesons melt [24| . The holographic description of the 
melting corresponds to the transition from the subcritical 
embedding to the supercritical one. 

Let us consider merger transitions in more detail. It is 
evident that the double-cone solution ([T]) is smooth ev- 
erywhere except for the tip p = where the curvature 
becomes singular. Near the tip the Kretchsman curva- 
ture scalar 1Z 2 = R^ vpa R^ pa infinitely grows 
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(2) 



The existence of the infinite curvature indicates that the 
solution obtained in the framework of classical Einstein 
gravity should be modified by quantum corrections. In 
other words, the naked singularity that is formed during 
the merger transition in its classical description might 
be resolved by the inclusion of quantum corrections into 
the classical action 25]. This conclusion is important. 
It means that if the transition between a black hole and 
black string phases occurs through the merger transition, 
one can expect the formation of a region with very high 
(up to the Planckian) curvature in a system characterized 
by macroscopic parameters (size of extra dimensions). 
An important question is how quantum gravity effects 
modify an adopted picture of classical merger transition. 

Trying to answer this question one inevitably meets 
two difficult problems. One is of technical origin, namely, 
how the near-critical solutions for the merger problem are 
modified by quantum gravity corrections, for example, 
by adding to the Einstcin-Hilbert action quadratic cur- 
vature corrections which arise in one loop computations. 
One can expect that the corresponding corrections be- 
come important when the curvature near the tip reaches 
the Planckian scale. The other more difficult problem is 
the following. At the corresponding Planckian scale the 
higher loop quantum gravity terms might also become 
important. If this happens, it would indicate that a com- 
plete solution of the problem requires the summation of 
all quantum loops or the use of a more fundamental the- 
ory of gravity, such as string theory. All of the above 
makes the problem very complicated for analysis. 

For this reason it is interesting to analyze a much sim- 
pler BBH model which has qualitatively the same behav- 
ior as merger transitions. Its critical solutions also have 
curvature singularity at the cone tip where it touches the 
horizon. One can expect that adding terms quadratic in 
the extrinsic curvature to the classical DNG Lagrangian, 
which are analogous to local one loop corrections in quan- 
tum gravity, may "cure" this "disease." Such curvature 
corrections naturally arise as a result of the stiffness ef- 
fect [26j . In the case of strings they were suggested by 
Polyakov [27| • We can think about such terms as cor- 
rections that come from the finite thickness of the brane 
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which is ignored in the DNG action. The DNG action 
can be considered as the zeroth order in the expansion 
in the width over a typical length in the system [H, H^j • 
Usually the small parameter in such an expansion is the 
ratio of the thickness of the brane to the characteristic 
radius of the brane bending. It is instructive first to 
study effects connected with the leading order terms in 
this expansion. In this analogy the thickness of the brane 
plays the role similar to the Planck scale. Moreover, if 
one describes the brane as a special topologically stable 
solution of some nonlinear field theory, one may, in prin- 
ciple, answer not only the question of how the quadratic 
curvature corrections modify the near-critical solutions, 
but also investigate the complete field theoretical object 
behavior in the near-critical regime. Effectively this cor- 
responds to the summation of all the stiffness corrections. 

In the present paper we focus on the first problem, 
namely we will analyze how the lowest order stiffness 
corrections modify the phase transition in the BBH sys- 
tems. 

A good analogy that helps to understand the effect of 
stiffness terms is a stiff bar. Consider the bending of a 
stiff bar. When we take into account the effect of the 
stiffness of the bar, its bending costs energy. Hence we 
cannot bend the bar as much as we want and it would 
eventually break long before a sharp tip is created. The 
sharp tip corresponds to the singular critical solution of 
the DNG action. One might expect that inclusion of 
higher derivative terms to the Lagrangian would prevent 
the creation of such a singular solution and will form a 
first order phase transition long before. Indeed, as we 
will see this is what happens in the BBH system for the 
subcritical configuration. 

The higher derivative corrections to the BBH system 
can serve as a toy model for the singularity resolution 
of "small BHs." Small BHs are singular limits of BH 
parameters in which the horizon becomes singular (see 
for example [3(1 [3l|). If we look at the induced BH on 
the brane, the critical solution has a singularity exactly 
of this type. 

The paper is organized as follows. In Sec. 77 we re- 
view the main results concerning the near-critical branes 
obtained in the absence of stiffness in the BBH model. 
In Sec. Ill we discuss the curvature corrections for a stiff 
brane. The stiff brane equations are presented in Sec. IV. 
Section V contains the analysis of near-critical solutions 
in the linear approximation. In Sees. VI and VII the 
numerical results for near-critical branes are presented. 
Section VIII contains a summary of the obtained results 
and their discussion. 



II. NONSTIFF BRANES 

In this section we briefly review the main results of [HI ] 
concerning the behavior of near-critical JJ-dimensional 
branes without stiffness interacting with a bulk static 
spherically symmetric TV-dimensional black hole. We do 



this mainly to explain the set up of the problem and to 
fix the notations we will use later. The Schwarzschild- 
Tangherlini metric of the bulk TV-dimensional spacetime 



dS 2 = g^daPda? = —FdT 2 + F -1 dr a 



2 dn 



N—2 



, (3) 



where F = 1 — (r g /r) 3 and d£l 2 N2 is the metric of 
a (TV — 2)-dimensional unit sphere S N ~ 2 . We define the 
coordinates 0j (i — 1, . . . , TV — 2) on this sphere by the 
relations 



dQ 2 +1 = d6 2 +1 + sin 2 e i+ idfl 2 



(4) 



We denote by x M (/i = 0, . . . , TV — 1) the bulk spacetime 
coordinates and by C a (a = 0, ...,£)— 1) the coordinates 
on the brane worldsheet. The functions x M = X M (£ a ) 
determine the brane world sheet describing the embed- 
ding of the (D — 1 )-dimensional object (brane) in a bulk 
TV-dimensional spacetime. We assume that D < TV — 1. 

In the absence of stiffness, the brane configuration in 
an external gravitational field g M „ can be obtained by 
solving the equations which follow from the DNG ac- 
tion [al, m H|| 



3 = / d D ^-&ct lab . 



(•5) 



where 7„b = g^X^X 1 ^ is an induced metric on the brane 
world sheet. We set the brane tension factor, which does 
not enter the brane equations, equal to 1. It is well known 
that an extremum of this action is a minimal surface. Let 



W 



be unit normals to the brane, and 



K 



(i) 
aj3 



(6) 



be an extrinsic curvature tensor. (V„ is a covariant 
derivative with respect to the bulk metric g^v) Then 
the nonstiff brane equations are of the form 



K {i) =g <*P K f Q . 



(7) 



For the axially symmetric D-dimensional static brane 
(with the isometry group 0(D — 1)) the induced metric 
is (n = D - 2) 



ds z = ~f ab dCd( b 

-FdT 2 + [F- 1 + r 2 (d9/dr) 2 ]dr 2 + r 2 sin 2 0dfi* , 
and the action ([5]) reduces to 

S = ATAn J drC, C = r n sin" 0y/l + Fr 2 {d9/dr) 2 . 

(9) 

Here AT is the interval of time, and A n = 2W 2 /T(n/2) 
is the surface area of a unit n-dimensional sphere. 



(8) 
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By analyzing the brane equation it is easy to show [16[ 
that for a brane which asymptotically approaches the 
equatorial plane 9 — n/2 one has 



p 

? = - 

r 



In? 



for n = 1 
for n > 1 



(10) 

We call the set of parameters {p,p'}, which characterizes 
the solution, the asymptotic data. 

The same solution can be determined by its behavior 
near the horizon. A subcritical brane is uniquely spec- 
ified by the distance of its tip from the horizon. The 
condition of the brane surface regularity at this point re- 
quires that its tangent plane at the tip is orthogonal to 
the symmetry axis. This fixes the second constant in the 
solution. Similarly, a regular brane crossing the horizon 
is orthogonal to the horizon surface, so that a unique 
constant fixing the solution is the "gravitational" radius 
of the induced BBH. A solution separating sub- and su- 
percritical solution is a critical solution. We denote by 
{p*iP'*} its asymptotic data. 



Asymptotic 
region {p,p'} 




Near (Rindler) 
zone (scaling 
transformations 
are valid) 



FIG. 2: This figure schematically shows a configuration of a 
supercritical brane in the regime when it is close to the critical 
one. Z is a proper distance from the horizon as a function of 
the radius R. Ro is the radius of the surface of the intersection 
of the brane with the horizon. In the region where Z <H r g the 
curvature surface of the horizon can be neglected (the Rindler 
domain) . 



Figure [2] illustrates the near-critical behavior of a su- 
percritical brane. (A similar graph for subcritical branes 
can be easily obtained from this one by evident changes.) 
A near-critical brane configuration is characterized by a 
parameter Rq, which is its radius at the intersection with 
the horizon. For a subcritical brane a similar parameter 
is Zq, the proper distance of the tip of the brane from 
the horizon. We consider a case when Rq (Zq) is much 
smaller than the gravitational radius r g of the bulk black 
hole. In the vicinity of the horizon, located at Z = 0, 
that is for Z<r 9l one has 



kZ 2 /2, F: 



k 2 Z 2 



(11) 



where n = ^(dF/dr)\ rg is the surface gravity. We call 
this region a near ( or Rindler) zone. The corresponding 
induced metric for a near-critical brane in the Rindler 
zone is 



ds 2 = -K 2 Z 2 dT 2 



dZ 
~dX 



dR 
~dX 



dX 2 + R 2 dn 2 n . 

(12) 

Here (Z(X), R(X)) is a brane equation written in a para- 
metric form. The action ((5|) for this induced metric is 



S 
S 



KATA n S , 

dXZR n ^(dZ/dX) 2 + {dR/dXf 



(13) 
(14) 



This action is evidently invariant under the transforma- 
tions A — > A(A). In the regions where either Z or R is a 
monotonic function of A, these functions themselves can 
be used as parameters. As a result, one obtains two other 
forms of the action which are equivalent to S 



S 



dZC, 



dRC 2 



(15) 



where 

C R = ZR n Vl + R' 2 , £ z = ZR n Vl + Z 2 . (16) 

Here the prime stands for the derivative with respect to 
Z, while the dot stands for the derivative with respect to 
R. The corresponding Euler-Lagrange equations are 

ZRR" + (RR' - nZ)(l + R' 2 ) = , (17) 



RZZ + (nZZ - R)(l + Z 



y2\ 



0. 



(18) 



It is easy to check that the form of Eqs. (fT"7|) - (fl"8|) is 
invariant under the following transformations: 

R{Z) = kR(Z), Z = kZ, (19) 
Z(R) = kZ(R), R — kR. (20) 

Equations ([TT}1 - ([T5 ]) have a simple solution 

R = VnZ (21) 

which plays a special role. We call it a critical solution. 
It describes a critical brane which touches the horizon of 
the bulk black hole at one point, Z = R = 0. It separates 
the two different families of solutions, supercritical and 
subcritical. 

The relation between Rq and is of the form 



In R = 7 In Ap + f (In Ap) 



(22) 



where 
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Ap = v/(p-^) 2 + (p'-p'J 2 , (23) 



and the function f(z) is periodic, f(z + lo) = f(z), with 
the period 



?r(n + 2) 
\[i + 4n - i 



(24) 
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Our aim is to study how the near-critical solutions are 
modified when a brane is stiff. We assume that the effec- 
tive width of the brane is much smaller that the gravita- 
tional radius of the black hole. In this approximation, as 
in the case of a nonstiff brane, the main features of the 
phase transition in the BBH system are determined by 
the brane behavior in the near zone, that is close to the 
event horizon of the bulk black hole, where the Rindlcr 
approximation is valid. 



III. STIFF BRANES 

In this work we consider the Dirac-Nambu-Goto action 
with minimal stiffness correction terms which play the 
role of higher curvature corrections [35| : 



jn+2 



CV^(l + B K 2 + CK 2 ). 



(25) 



Here 7 is the determinant of the induced metric given in 
Eq. ([5]), K = J2i Kjf 1 ^ is the trace of the extrinsic cur- 
vature tensor, and K, 2 — J2i K(i)ni/K " is its square. 
A minimal model of stiffness corrections involves only 
quadratic powers of the extrinsic curvature tensor. In 
this work we will concentrate on the above "truncated" 
model as a toy model whose solution gives us the static 
configuration of a stiff brane embedded close to the hori- 
zon of the bulk black hole, namely, in the Rindler zone. 

In the particular case of a domain wall the stiffness 
coefficients were calculated in the framework of a mi- 
croscopic model of a vacuum defect in field theory with 
spontaneous symmetry breaking (29| . In this work we 
show that the exact numerical values of the coefficients 
do not affect the qualitative features of the solution. Nev- 
ertheless the sign of the coefficients is important. Positive 
stiffness coefficients 

B,C > 

ensure us that the energy density for the static solution 



e = -£ = v/=7(l + BK 2 + CK?), 



(26) 



is positive since K 2 and K? are both non-negative. 

Before going further let us discuss two interesting spe- 
cial cases of the general theory. 

C = case. In this special case the stiff string equa- 
tions have a simple exact solution. Namely, any solution 
of the DNG equations ([7]) is at the same time a solution 
of the stiff string equations. Indeed, a general variation of 
the action can be split into the variation along the brane 
and the transverse one. Variations along the brane sur- 
face vanish identically. For transverse variations, multi- 
plying the equations of motion by a normal to the brane 
gives 



-2BJ 



-7 K n 



(l + BK 2 ) 
SK 



= 0. 



(27) 



Now, substituting the DNG equation K^> = into the 
right hand-side of (f27|) . we see that it vanishes. Hence 
i£w — is a solution of the stiff string equations. 

B + C = case. This case is not interesting for our 
consideration since it violates the positive energy condi- 
tion, but it is of mathematical interest. Let us notice 
that the Gauss-Codazzi relation for a flat bulk spacetime 
implies 



K = K 1 - K 2 



(28) 



Thus one can rewrite the action (|25|) in the following 
form: 



jn+2, 



BK+{B + C) K? 



(29) 



For B + C = the term with JC 2 vanishes. 

Let us return to the discussion of the general case. 
The units of the stiffness coefficients are length-squared. 
Therefore under scaling transformation of the spatial co- 
ordinates 



R — ► s R , Z ^ sZ, 



we have 



'2 , „ — 2 r/2 i^2 , s ~ 2 jq2 



K 2 -> s- 2 K 2 , K 2 



(30) 



(31) 



Using this transformation we can set one of the coef- 
ficients to be unit, say C = 1. We can think about it 
as taking the basic unit length of the stiff string to be 
y/C. Then we are left only with one free parameter B in 
the action. We will use this choice later in Sees. I VII and 
IVII1 when we will discuss the results of the numerical 



calculations. 

For completeness we give here the components of the 
extrinsic curvature for the brane (Z(X),R(X)) with the 
induced metric (TTJ 



Ktt 



V~ X A, 
dR 

- z dx v ~ 

dX ' 



(32) 



where 



V 
A 



xJ{dZ/dX) 2 + (dR/dX) 2 
dZ d 2 R dR d 2 Z 

~dx ~d>? ~ ~d~x Tlx?' 



(33) 
(34) 



The action (f2"5| for the induced metric (fT2")) takes the 
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following form: 



s 
c 

Co 

c 2 



-kATA„ / dXC 



(35) 



Co + Bd +CC 2 
ZR n V, 



ZR-V 

R n {dR/dX) 2 

zr 

ZR n A 2 



dR/dX 



n dZ/dX \ ' 



ZV RV 

nZ{dZ/dX) 2 R n - 2 



•p5 



The action (|25|) is evidently invariant under transforma- 
tions A — > A(A) as in the nonstiff case. In a general case, a 
variation of this action gives equations containing fourth 
derivatives, while the corresponding constraint equations 
are of the third order I39| . 



IV. EULER-LAGRANGE EQUATIONS FOR 
STIFF BRANES 



In the regions where either Z or R is a monotonic func- 
tion of A, one of the coordinates can be used as a param- 
eter. As a result, one obtains two additional forms of the 
action: 



S = - / dZCn = - / dRC 



where 



C R = ZR n V \l+B 



-C 



112 



R" R' : n 
V s + ZV + ~RV 

R' 2 n 



R 



v e Z 2 V 2 



R 2 V 2 



C z = ZR n V \1 + B[ 



7 2 



ZV 



Z 
pa 



Z 2 V' 2 



n Z 
~RT^ 

nZ 2 > 

rJv 2 ^ 



(36) 



(37) 



(38) 



V in the first relation means V + i?' 2 , while in the 

second equation one has 

As in the nonstiff case, the form with R(Z) is useful 
for the description of the supercritical brane while Z(R) 
is more suitable for the description of subcritical branes. 

The equation for R(Z) takes the following form: 



- 2 b Z 3 R 3 (1 + R' 2 ) 2 i? (4) + 4 b Z 2 R 2 (R' 2 + l) [5 Z RR' R" - (R' 2 + l) (n Z R' + R)] R {3) + F (R, R , R" , Z) = 0, 

(39) 

where 



F(R,R',R",Z) = bbZ 3 R 3 (l -6R' 2 ) R" 3 + 3bZ 2 R 2 (l + R' 2 ) [5RR' + nZ (AR 12 - 1)] R" 2 
Z 3 R 3 (l + i?' 2 ) 3 R" -ZR (l + R 2 f (2 [2b + 3B] nZ RR' + bR 2 [R 2 - 2] 

nZ 2 [b + 3 B - 3 B n + 2b (n - 2) R 2 ]) R" + (l + R' 2 f (R 3 R [Z 2 + Z 2 R' 2 ] -nZR 2 [Z 2 + Z 2 R' 2 ]) 



(40) 



- (1 + R' 2 ) 3 ([b-3B(n- 1)] nZ 2 RR' - bn Z R 2 R' 2 + b R 3 R' [R 12 + 2] +n[n-2] Z 3 [b + B (n - 1) + 2 bR' 2 ] ) , 

and b = B + C, 

A similar equation for Z(R) is 

-2bR 3 Z 3 (l + Z 2 Y Z (i) +4bR 2 Z 2 (l + Z 2 ^ 5RZZZ- (l + Z 2 ^ (nZ + RZj Z (3) + g(z,Z,Z,R) = 0, 

(41) 



where 

g(z,Z,Z,r) = hbR 3 Z 3 (1- 6-2*1 Z A + 3bR 2 Z 2 ll + Z'" 



r,nZ Z + R (aZ 2 -1 



Z 2 + R 3 Z 3 ( 1 + Z 2 ) Z 



72 



RZ[\ + Z 2 ) (2[3B + 2b]nRZZ-bR 2 



2Z 2 -I 



nZ 2 [2b{n-2) + ('3B + b-3B ri)] ) Z 



R 2 Z 2 (R-nZZj (l + Z 2 ^j A + (l + Z 2 Y (^-bnR 2 Z Z + nRZ 2 [6 - 3B(n - 1)] Z 2 + bR 3 



1 + 2Z 2 



(42) 



2 b n (n - 2) Z 3 Z + n (n - 2) Z 3 Z 3 [b + B{n - 1)] 



Let us denote 



I has dimensionality of the length. The extrinsic cur- 



/ = max ( s/B, VC 



(43) 
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vature corrections are dominant for R, Z < I where the 
stiff branc differs significantly from the DNG brane. The 
significant effect of the stiffness is localized in the region 
where the original DNG brane is extremely bent. This 
happens in the neighborhood of the point R = Z = 
which is defined by the length scale For R, Z ^> I 
the solution for the stiff brane-BH system approaches 
the DNG brane-BH system. In particular, R = y/nZ 
is the attractor solution for the DNG brane-BH system 
and therefore it should be also an attractor for the stiff 
brane-BH system. 



V. NEAR-CRITICAL MODES 

Let us study linear perturbations to the attractor so- 
lution for the case of a stiff brane. Our objective is to 
obtain the modes in the neighborhood of the attractor 
far away from the singular region Z 3> I but still located 
in the Rindler zone. The modes will guide us later in the 
setting of the boundary conditions for stiff branes. 

Let us substitute in the equation for R the follow- 
ing expression: 



R{Z) = y/nZ + p{Z), 



(44) 



and keep only linear terms in p(Z). Then we obtain the 
following linearized equation: 

a%{Z)p + a Q l s (Z)p' + a s 2 {Z)p" + a s 3 (Z)pW+al(Z)pW 
= 2C(n- I) (n + ifrT^Z, (45) 

a s (Z) = (n + 1) [25 (n- 5) - C(n + 7)] + a (Z), 

a{{Z) = (n+ l)Z[2B (16 -n) - 17C(n- I)] +a 1 {Z), 

a s 2 (Z) = -(n + l)Z 2 [2B(n + l) + C(2n-l)] + a 2 (Z) ) 

al{Z) = -A{B + C){n + l)Z 3 , 

a\{Z) = -2{B + C)Z A . (46) 

Here a,i(Z) are the coefficients in the linearized DNG 
brane equations: 



a Q (Z) = {n+l) 2 Z 2 , 
Oi(Z) = {n+l) 2 Z 3 , 
a 2 {Z) = (n+l)Z 4 . 



(47) 



Now let us take the limit Z ^> I and as a result we obtain 
the following equation: 

(n + 1) 2 [Zp + Z 2 p 1 ] + (n + l)Z 3 p" 

- 2{B + C) |2(n + 1) Z 2 p {3) + Z 3 p (4) 

= 2C{n- 1) (n+ l) 2 n"5. (48) 

The leading term of the particular solution at large Z 



is 



PP 



C{n 2 



hiZ 



O 



1 

Z* 



and therefore it does not have an effect on the attractor 
R = \/nZ at large Z, as expected. A general solution 
of (|4"B")) is a sum of this particular solution and a general 
solution of the homogeneous equation obtained from (f4"5|) 
by omitting the right-hand side. It is interesting that this 
homogeneous solution depends only on the sum B + C of 
the stiffness coefficients. 

Since the homogeneous equation is of the fourth or- 
der, it has four linearly independent solutions (asymp- 
totic modes). Two of the asymptotic modes reproduce 
the asymptotic solutions for the DNG brane: 



p ^ Z -i(n±Vn*-4n-4)_ ^ 

The other two modes appear only for the stiff brane: 



exp ±, 



1 



2 (B + C) 



Z 



(51) 



(49) 



The additional two modes above are added due to the 
stiffness corrections. One of the additional modes is an 
unstable mode which takes the solution away from the 
attractor. This mode should be eliminated by appropri- 
ate boundary conditions in order to reproduce the DNG 
solutions at large distances. 

Thus we arrive to a boundary value problem. Let us 
take, for example, the subcritical configuration where the 
solution can be written as Z{R) (a similar discussion is 
applicable to the supercritical configuration with some 
evident changes). Since Eq. (|4"Tj) is of the fourth order, 
we need four initial values (in case of an initial value 
problem). The configuration is axially symmetric with 
the symmetry axis R = and it is plausible that the 
stiff brane solution preserves the same axial symmetry. 
Consider a brane passing through the point Z(Q) = Z 
(the proper distance of the tip of the brane from the 
horizon). The axial symmetry and the regularity of the 
brane at R = enforces Z(0) = and Z^{0) = 0, while 
Z(0) remains a free parameter. This free parameter will 
allow us to eliminate the unstable mode by finding the 
specific value for Z(0). 

In conclusion, we have a boundary value problem defin- 
ing near-critical solutions of the stiff brane equations: In 
order to obtain the right asymptotic behavior at large R 

R 



for a brane which passes at Z(0) — Zq we have to tune 
the parameter Z(0) = Zq. 



VI. NUMERICAL RESULTS FOR STIFF 
BRANES: n = 1 CASE 

Using a numerical shooting analysis we obtained the 
values of Zq for which there is a solution for the sub- 
critical brane that satisfies the boundary conditions: It 
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starts at Zq and asymptotically goes to the attractor. A 
similar analysis was performed for the supercritical con- 
figuration where values of R"(Q) were determined as a 
function of Rq (the radius of the BBH horizon) . 

Let us start by examining the case n = 1. As we will 
see, this case is qualitatively different from n > 1. For 
B = the action [Eq. ([35]) ] is completely symmetric for 
the interchange of R <-> Z. This discrete symmetry of 
the equations implies the same results for the subcritical 
and supercritical configurations. For this reason we give 
here the results for the subcritical configuration B = 
and n = 1 (Fig. [3]) when for the supercritical configu- 
ration the graph is the same (up to the interchange of 
R «-> Z). The plot at this figure shows Zq as a function 
of Zq for near-critical configurations. The finite gap in 
the neighborhood of the point R = Z = demonstrates 
the first order phase transition in the BBH system. A 
detailed interpretation of this picture is given below in 
the discussion of the case n = 2 (where the same picture 
appears only in the subcritical configuration). 



Zo 

0.08 
0.06 
0.04 
0.02 



B = 



J 2 3 4 5 6 



Zo 



FIG. 3: as a function of Zo for n = 1 and 5 = 0. The 
symmetry of the action in this case implies that the same 
graph is valid to the supercritical configuration as well - R" (0) 
as a function of Rq. 



into Eq. (|4"Tj) and keeping only linear terms in C(R) gives 
the following linear equation (n — 1) 



(B + C) i? 4 C > (R) + 4 (B + C) R A C > (R) 
R 2 (AB + C- R 2 ) C(R) -2R 3 ((R) 
2(R 2 -AB-AC) C(R) =0. 



(55) 



Hence the symmetry R <-> Z is demonstrated analyti- 
cally in the linear approximation. This does not imply an 
exact reflection symmetry of the solutions, but at least 
makes it possible. 



VII. NUMERICAL RESULTS FOR STIFF 
BRANES: n > 1 CASE 

For n > 1 there is no reflection symmetry of the ac- 
tion anymore, and sub- and supercritical solutions be- 
have quite differently. Let us start with the subcritical 
configuration and demonstrate that it exhibits the same 
qualitative features as n = 1. This is a good place to 
compare the details of the new picture with the nonstiff 
case. 

Consider, for example, the case of n = 2 and B = 1 
in Fig. |4] The value of Zq is plotted as a function of the 
position where the brane crosses the axis of symmetry 
Zq. In the case of DNG branes, i.e. without stiffness, the 
dependence of Zq on Zq is determined from the Euler- 
Lagrange equation (fT8| to be (see in [161]): 



Zq 



1 



(n + 1) Zq 



(56) 



This function is plotted in Fig. 
the case of stiff branes. 



[4] for comparison with 



Zo 



For the case of B ^ the action (|35j) is no longer sym- 
metric under the reflection R <-> Z. Nevertheless we find 
numerically that the results are symmetric within the 
used accuracy. Evidence for this symmetry we can find 
in the linearized equations for a perturbation around the 
attractor. In Sec. |V]we studied linearized perturbations 
to the supercritical configuration 



R{Z) = ^Z + p{Z). 



(52) 



Keeping only linear terms in p(Z) gave us Eq. (|45[) . For 
the special case of n = 1 this equation reads 

(B + C) Z 4 (Z) + 4(B + C) Z 3 p ( V (Z) 
+ Z 2 (AB + C - Z 2 ) p"{Z) - 2Z 3 p'(Z) 
- 2 (Z 2 -AB-AC) p{Z) = (53) 

In a similar way for the subcritical configuration sub- 
stitution of 



Z(R) = ^=+t(R) 



(54) 



B = 1 



0.15 



0.1 



0.05 



Zo 



FIG. 4: Zq as a function of Zq for n = 2. The dashed line is 
the same function for DNG branes (without stiffness terms). 

Few features can be observed in the graph: 

• There is a finite gap < Zq < 1 in which the 
solution for the embedded stiff brane does not exist 
at all. Hence the singular point is resolved for the 
subcritical branch. This is a characteristic feature 
of first order phase transitions. 



• Zq is bounded, unlike DNG branes [Eq. 
which Zq is unbounded 



for 



Zo 



B = 100 



For 1 < Z < 1.25 we see coexistence of two 
branches of solutions. For any Zq in this range 
there are two possible values of Zq and thus two 
possible configurations of the stiff brane. One can 
compare the energy [see (1551) ] of the two branches. 
A numerical comparison of the energies [Fig. [5] re- 
veals that the branch with the lower values of Zq 
is energetically favored. The branch with higher 
energy corresponds to a local phase at maximum. 
This solution should be unstable and separates two 
stable phases in a first order phase transition. 



205 
20(1 
195 
190 
185 
180 



Zo 



FIG. 5: The energy density integrated for < R < 5 as a 
function of Zo comparing two branches in the segment (1 < 
Zo < 1.25). Note that the minimal energy is obtained at the 
point which corresponds approximately to Zo =0. 
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FIG. 6: Z as a function of Z for n = 2. B = 100, 0.01. 
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0.1 
0.05 



B = 1 



Zo 



• There are solutions for stiff branes that satisfy the 
boundary conditions with negative Zq. Such solu- 
tions exist only for -0.025 < Zq < 0. 

• There exists an "end point" in the plot with min- 
imal value of Zq. For Zq less than this value a 
solution does not exist. 

• At large values of Zq we see that the effects of stiff- 
ness are negligible. The points of the stiff branes 
approach the DNG branes at large values of Zq. 

In order to check that the obtained results are robust 
we repeated the same calculations for various values of B. 
In all cases we found that the same qualitative behavior 
repeats itself: A finite gap in the existence of solutions 
for < Zq < I, two branches of solutions in a small 
neighborhood of Zq ~ / etc. 

For illustration we give in Fig. [6] two graphs for two 
values of B with four orders of magnitude difference (B = 
0.01, 100). 

In addition we checked for various dimensions n = 4, 5 
and found the same qualitative behavior (see Fig. [7] for 
n = 4 as an example). Despite the fact that n < 4 
is different from n > 5 since in the former the phase 
space behavior of the critical solution behaves as of a 
focal point and in the latter as a node (see [IB]). This 
type of transition in the near-critical solutions has no 



FIG. 7: Zo as a function of Zo for n = 4, B = 1. 



influence on the neighborhood of the singular point R = 
Z = where the stiffness terms are dominant. 

It is surprising that when we repeated similar calcula- 
tions for supercritical stiff branes with n > 1 we found 
that for the supercritical configurations there is no singu- 
larity resolution. The stiffness terms break the symme- 
try between the supercritical and subcritical brane-black 
hole systems. The supercritical solutions show no gap nor 
double-branch behavior. As an example let us look at the 
supercritical configuration for n = 2. For B < 0.906 we 
did not find evidence for the existence of a solution in 
the vicinity of the point R(Q) = R"(0) = 0. See Fig. M 
for B = 1 and Fig. M for B = 0.5. We stress that in both 
cases the curvature singularity still exists. 



VIII. DISCUSSION 

We observed that due to the stiffness corrections the 
singularity of the critical solution is resolved for n = 1 
in a symmetric form (both in the subcritical and super- 
critical configurations) and for n > 1 only for subcritical 
configurations. We observed this resolution in the cre- 
ation of a finite gap and a clear signature of a first order 
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FIG. 8: R"(0) as a function of Rq (supercritical) for n = 
2, 5 = 1. 
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FIG. 9: R"(0) as a function of Ro (supercritical) for n = 
2, B = 0.5. 



merger transitions when higher derivative corrections are 
included. Inclusion of higher derivative corrections might 
cause the merger transition to become first order in na- 
ture and create a finite gap between the thin black string 
("the waist") and the caged black hole. This way the 
naked singularity and the violation of cosmic censorship 
hypothesis that appear in the classical approximation 
would be resolved. This might also be a natural way to 
resolve the apparent tension between the suggested sce- 
nario for the merger transition and the observation that 
such a pinch-off can occur only at infinite affine param- 
eter along the horizon [3?|]. The resolution is the follow- 
ing. When the system approaches the Planckian scale, 
at a finite time, the first order phase transition takes the 
system to the second phase. Therefore with quantum 
corrections the "pathologies" of infinite affine parameter 
and naked singularity would be resolved. 

The asymmetry that we found might be a result of the 
incompleteness of the truncated model that we used to 
describe the full effect of quantum corrections. It might 
be also a hint on an asymmetry which is generic in the 
topology change in general. 
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